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In order to control aggregation phenomena in suspensions of nanoparticles, one often charges the particles electri- 
cally, e.g. by triboelectric charging. Stabilization of suspensions against aggregation of particles is an important 
issue, which may be realized by monopolar charging, where particles repel each other. Contrarily, bipolar charging 
may be used in coating processes, where smaller particles of one material coat larger particles of another material. 
When the two particle fractions are charged oppositely, aggregation between equally charged particles is hindered 
whereas aggregation between oppositely charged particles is preferred, thereby improving the coating process. 
We study various aspects of these two situations by theoretical investigations and computer simulations. 



1 INTRODUCTION 

Granular materials play a major role in a large num- 
ber of technological applications. The size of grains 
used depends on the application, ranging from large 
particles with diameters of several centimeters down 
to nanoparticles. However, the physical properties of 
powders change drastically when reaching the size of 
nanoparticles. While in ordinary dry granular materi- 
als contacts between the particles cannot sustain ten- 
sile loads, in the case of ultra-fine powders van der 
Waals forces lead to aggregation which is essentially 
irreversible. In contrast to volume forces such as grav- 
ity, which decrease with decreasing particle radius a 
like a 3 , van der Waals forces are proportional to a 
and thus dominate all other forces in the case of suf- 
ficiently small particles. On the other hand, van der 
Waals forces are only short ranged, diverging like h~ 2 
on very short distances and decaying like hr 1 with 
increasin g distance h between the surface of parti- 
cles ( IT. van de Ven 19891 . Thus, the van der Waals 
force can be approximated as an irreversible sticking 
force, which has no influence on particles separated 
from each other, but causes them to stick together ir- 
reversibly, as soon as they touch each other. 

A promising strategy to reduce agglomeration of 
particles is to charge them electrically. If one man- 
ages to charge all particles with the same polarity, the 
mutual Coulomb repulsion delays or even suppresses 
the aggregation process. It has been shown that such a 



charge distribution can be establis hed for particles sus- 
pende d in a non- polar liqui d O.H. Werth et al. 20031 
M. Linsenbiihler et al. 2005 1. In the first part of this pa- 
per we investigate to what extent charges can be used 
to enhance the stability of a suspension against ag- 
glomeration. We present some results obtained from 
a rate equation investigation, rec ently published in 
iSM. Dammer & D.E. Wolf 20041 

In the second part of the paper, we present an appli- 
cation which makes use of Coulomb forces: the coat- 
ing of charged particles with sizes on the micrometer 
scale by oppositely charged nanoparticl es. As shown in 
iM. L insenbiih ler & K.-E. Wi rth 2004 1, such a coating 
may enhance the flowability of a powder dramatically, 
since the nanoparticles on the particle surfaces act as 
spacers between the micrometer sized particles, thus 
reducing their mutual van der Waals forces. The aggre- 
gation process of oppositely charged particles has been 
studied in detail (U.H. Werth et al 20 05 1. 

2 STABILIZATION AGAINST AGGLOMERA- 
TION 

Let us consider a very fine suspended powder. The 
size of the particles shall be a few /xm or below. We 
assume the particles to be all of the same insulating 
material, each carrying an electric charge of the same 
sign. The charged particles are subjected to Coulomb 
and van der Waals forces as well as Brownian mo- 
tion and Stokes friction. As already pointed out in 



the introduction, van der Waals forces will only dom- 
inate over Coulomb forces for particles very close to 
each other. Thus we assume that the van der Waals 
force has no effect except when particles collide and 
stick together irreversibly. The motion of the particles 
is then composed of a stochastic part due to Brow- 
nian motion and a deterministic part determined by 
the balan ce of repelling Cou lomb forces and Stokes 
friction ( S.M. Dam mer et al. 2 004 1. It may happen that 
particles (or agglomerates of particles) come into con- 
tact due to Brownian motion, thus forming a new ag- 
glomerate (Brownian coagulation). In the following we 
will investigate this aggregation process in further de- 
tail in terms of a rate equation approach. 

We assume that initially the suspension consists of 
primary particles or primary aggregates, all of the same 
mass m*, the same effective radius a*, and carrying the 
same charge q* . As particles collide they form clusters 
with increasing mass and charge. Since for the chosen 
initial conditions mass and charge of the clusters are 
proportional to each other, it is sufficient to describe 
them by a single index i according to 
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Here m (t) denotes the number density of clusters with 
mass rrii—im* and radius ai—i a a* at time t, each of 
them carrying the charge qi=iq*, where 1/a denotes 
the fractal dimension of the aggregates (e.g. a=l/3 
for spherical particles). As initial condition we choose 
ni(t=0)=l, Ui (t=0)=0 for i>l . The matrix Rij in 
the rate equation is called coagulation kernel and de- 
scribes at which rate two clusters with indices i and j 
merge into a single one. For Brownian coagulation of 
charged particles one has 



k 2 ij 



exp(Ky) - 1 



(i a +j a ) 



,(2) 



(3) 



where k 2 =q* 2 /{Altera* h^T) and 2fc B T/3r/=l 
(IS.M. Dammer & D.E. Wolf 20041 . Here n t3 is pro- 
portional to the Coulomb energy of two clusters 
being in contact, divided by the thermal energy. This 
Coulomb energy is necessary to bring two particles 
into contact from infinite separation. If the suspension 
contains countercharges so that a double layer forms, 
the Coulomb repulsion at short distances is combined 
with attractive forces at intermediate distances. In this 
case the activation energy for a particle collision can 
be enhanced, i.e. the collision rate can be reduced 
compared to Eq. (0. 




Figure 1. The average mass M(t) scaled with M c vs. time 
t scaled with t c for different values of k. The data collapse 
shows that in all cases the crossover to slow cluster growth 
happens at t=t c and M=M C . The inset shows the original 
data, and the dashed line is the behavior for an uncharged 
system (k—0). 

If the system is initially unstable, i.e. fc 2 <Cl, 
it behaves essentially as if it was uncharged and 
aggregation occurs frequently. In this regime, 
the average aggregate mass increases linearly 
in time JM.H. Ernst 19861 . However, as soon 
as particles with Kijfxl become important, fur- 
ther aggregation is suppressed exponentially. 
Asymptotically, in this regime the average ag- 
gregate ma ss increases sub-logarithmically slow 
in time (IS.M. Dammer & D.E7Wolf 20041 The 
crossover between these two regimes happens at a 
characteristic time t c and mass M c given by 
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where t* is the appropriate time unit, given by t* = 
3r]/h2ksT. h is the initial number density of pri- 
mary particles. Fig. \l\ shows the temporal evolution 
of the average aggregate mass, obtained from a nu- 
merical solution of the rate equations Q. As an ex- 
ample, for an experimental situations with values like 
J7=10 _4 Pas, T=300K, a*=l^m and n«2-10 5 mm- 3 
one finds £*«150s. Assuming e r «l, charging with a 
single elementary charge on each primary particle cor- 
responds to a value fc«0.2. For these parameters one 
finds with the dimensionless results of Fig. [2 that the 
crossover to slow aggregation occurs within hours. 

For times larger than t c the suspension is usu- 
ally called stable. However, note that aggregation still 
proceeds slowly. A consequence of this is that the 
mass distribution evolves towards a universal scal- 
ing form and that the relative width ay of the dis- 
tribution falls off to a universal value <jj?°«0.2017 
that is much smaller than in the uncharge d case (with 
oywl) (IS.M. Dammer & D.E. Wolf 20041 Hence, for 
an initially narrow mass distribution ay first grows sim- 
ilarly to the uncharged situation until time t c , when 
it starts decreasing again (induced by the exponential 
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suppression of further aggregation) which is denoted 
as self-focusing. With the values of the example dis- 
cussed above, this self-focusing effect occurs within an 
experimentally accessible time. 

3 COATING OF PARTICLES WITH NANOPARTI- 
CLES 

In this section we present an application which 
makes use of charged granular particles in sus- 
pension, namely the coating of powder parti- 
cles with nanoparticles. As shown b y experiments 
iM. Linsenbiihler & K.-E. Wirth 2 004 1, such a coat- 
ing of powder particles, even by a small amount of 
nanoparticles, leads to a strong increase in flowability 
and fluidizability of powders. 

A detailed description of the coating process can 
be found in JJ.H. Werth et al. 2003> . Therefore, we 
want to give only a brief description. During the pro- 
cess, both the powder particles and the nanoparticles 
are charged triboelectrically by stirring the suspension 
with a high speed stirring device. By choosing dif- 
ferent, electrically insulating materials for the pow- 
der particles and the nanoparticles it is possible to 
charge both fractions oppositely in a single process. 
After the charging process is stopped, the suspension 
is left untouched and the particles may agglomerate 
with each other. Due to Coulomb forces, agglomera- 
tion of oppositely charged particles (i.e. powder parti- 
cles and nanoparticles) is preferred, while agglomera- 
tion of equally charged particles is suppressed. Eventu- 
ally, a dry powder of coated particles can be achieved 
by evaporation of the liquid nitrogen. 

An interesting question is, to what extent the charges 
influence the deposition of the nanoparticles on the 
larger powder grains. To this end, we expect the 
nanoparticles to usually carry only one single elemen- 
tary charge, while the larger powder particles may 
carry several elementary charges. Since the used mate- 
rials are electric insulators, the elementary charges are 
immobilized at distinct positions on the particle sur- 
faces. If a nanoparticle now approaches the surface of 
a larger particle, two situations may arise: it may touch 
the particle exactly at the position of a charge, thereby 
compensating this charge, or it may just touch on a ran- 
dom position on the surface, leaving a dipole moment 
interacting with more charged nanoparticles in the sus- 
pension. The second case may occur due to Brownian 
motion of the particles. 

To answer the question, which situation is more 
likely and to what extent surface charges of powder 
particles are compensated by charged nanoparticles, 
we study a simplified model: point-like particles, each 
carrying one negative elementary charge, are conti- 
nously inserted far away from one infinitely large par- 
ticle, represented by a flat wall and carrying one or sev- 
eral positive elementary charges. The point-like par- 



ticles are subjected to Coulomb forces by interaction 
with the charges on the wall, to Brownian motion and 
to Stokes friction by the surrounding fluid. Coulomb 
interaction of equally charged particles among each 
other is not regarded. As soon as a particle touches the 
wall, it is removed from the system. Our aim is to com- 
pute the probability distribution p(r ) for the particles to 
touch the wall at position r. 

The answer to this problem is given in detail in 
JJ.H. Werth et al. 2*0 05 ). Therefore, here we only moti- 
vate the according equations and present some results. 
The probability distribution P(f, t) to find a particle at 
point r and time t can be described by a Fokker-Planck 
(FP) equation. It has the form 

^- t P(r,t) = -V-j(f,t), (5) 

where 

j (r, t) = -L>VP(r, t) + v{r)P(r, t) (6) 
is the probability current, D the diffusion constant, and 
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is the particle velocity in the overdamped limit. 

Measuring space in units of Q/D and time in units 
of Q 2 /D 3 we obtain the parameter-free dimensionless 
equation 

d -> 

— P = V 2 P-VP-u-P(V-u) (8) 
dt 

where u = —rjr i . 

We are only interested in the distribution p(r) of in- 
coming particles at the wall. This allows us to use a 
time independent FP equation. We arrive at the station- 
ary FP equation 

V 2 F- VP-u-P(V-u) = 0. (9) 

Equation can be solved analytically in the case 
of one single charge fixed at the wall. The analytical 
solution as well as an in-depth derivation of (|9} and a 
numerical treatment of the ac cording Langevin equa- 
tion can be found in dJ.H. Werth et al. 20 05 1. 

First we want to concentrate on the analytically solv- 
able case of one single surface charge at the wall. The 
density distribution p(r) of particles reaching the wall 
in a distance r to the fixed charge is given by 

p(r) = l + ±- + ^6(r), (10) 

where the density of nanoparticles far away from the 
wall is set to unity. There are three terms contribut- 
ing to the density distribution: Firstly, there is a cer- 
tain fraction of point particles which exactly reaches 
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line of charges by y, the light curves are given by 
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Figure 2. Density profile of particles hitting the wall. One 
(a), two (b), three (c) and four (d) charges are fixed on the 
wall in a line. The pictures show density profiles of incoming 
particles along this line. The bold curves present numerical 
results, while theoretical predictions are plotted in light grey. 

the surface charge, represented by the term §<5(r). Sec- 
ondly, there is a concentration of particles around the 
fixed charge, decaying like l/2r with growing distance 
r to the fixed charge. Finally, there is a constant back- 
ground of incoming particles according to the density 
of particles far away from the wall. Solution dlOt pro- 
vides no intrinsic length scale. 

Moving from one single fixed charge to several 
charges distributed at the wall, there is no analytical 
solution of equation (|9jl until now. Especially, the sum 
of solutions of type i ll Ob for several charges at different 
positions is not a solution of the multi-charge Fokker- 
Planck equation. However, as we will see, at least for 
only a few charges, adding up the single charge solu- 
tions can be a first approximation for the real problem. 

Since analytical solutions are not tangible, the di- 
rect numerical simulation of a multi-charge Langevin 
equation can be used to obtain results. The Langevin 
equation describing the movement of a nanoparticle is 
given by 
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where r is the position of the particle, Fc(r) is the 
Coulomb force acting on it, rj is the dynamic viscos- 
ity of the fluid, a the particle radius and £(t) is a white 
Gaussian noise corresponding to Brownian motion of 
the particle. 

Fig. |2] shows densities of particle influx at the wall 
obtained numerically. In the simulations corresponding 
to Figures |2] (a)-(d), up to four charges are arranged 
in a line, separated by a distance of 2Q/D. The bold 
curves show the density of particles touching the wall 
along this line. The light curves are obtained by adding 
up the l/2r-shoulders from the single charge solution 
(II Oi for all charges. If we denote the position along the 
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where n is the number of charges fixed at the wall. 
This approximation agrees nicely with simulations, al- 
though in the case of three and four charges the numer- 
ical data lie slightly below the light curves. 

4 CONCLUSION 

We discussed suspensions of monopolarly and 
bipolarly charged particles. In the former case, accu- 
mulation of charges leads to a crossover from fast to 
slow aggregation, followed by a self-focusing of the 
mass distribution towards a universal scaling form. 
In the latter case, we showed that bipolar charging 
may enhance the coating of powder particles with 
nanoparticles. 
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